Landau Theory of the Phase Transitions in Half Doped Manganites: Interplay of 

Magnetic, Charge and Structural Orders 
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The order parameters of the magnetic, charge and structural orders at half-doped manganites are 
identified. A corresponding Landau theory of the phase transitions is formulated. Many structural 
and thermodynamical behaviors are accounted for and clarified within the framework. In particular, 
the theory provides a unified picture for the scenario of the phase transitions and their nature with 
respect to the variation of the tolerance factor of the manganites. It also accounts for the origin of the 
incommensurate nature of the orbital order and its subsequently accompanying antiferromagnetic 
order. 
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The discovery of "colossal" magnetoresistance has 
stimulated aj-enaissance of interest in doped rare-earth 
manganese .EII3 Intensive investigation has revealed a di- 
versity of novel phenomena due to the complex interplay 
among magnetic, charge, orbital and structural orders. 
A particular relevant issue is the competition between 
magnetic and charge orders in half doped manganites. 
It poses a great challenge to theorists upon how to deal 
with the strong correlation in models with magnetic, or- 
bital and lattice degrees of freedom. Here we formulate a 
Landau theory of phase transitions based on the symme- 
try of the system in an attempt to understand a variety 
of sometime controversial structural and thermodynam- 
ical behaviors. Although it may be argued to be only a 
mean field theory which is incorrect at critical points, the 
structural information, among others it affords is robust. 
And it is the order parameters that exhibit singularity at 
the critical points. 

The most prominent charge ordered (CO) behavior in 
perovskite manganites concerns with those doped at 0.5. 
These systems exhibit several classes depending on the 
tolerance factor of the resultant structure (see Fig. |^). 
For Lao.sSro.sMnOa with small distortions, which we 
classify as Class I hereafter though no charge order ap- 
pears (neither does Class II below), a paramagnetic (PM) 
to ferromagnetic (FM) transition occurs at Tc ~ 360K. 
When La is replaced by a smaller ion Nd, Tc decreases 
with the tolerance factor. At 40 percent of Nd or so, an 
intermediately distorted class II shows up in which the 
FM phase transforms at a lower temperature Tajm/ to 
a metallic A-type antiferroMiagnetic (AFM) statdj^Q. As 
more La is replaced by Nd,El a new class III sets iruwhich 
displays a CO CE-type AFM state below Tco^ Sia^ 
ilar behavior has been reported in Lao.5Cao.5Mn03.El'El 
Pro.sSrg.sMnOs iSj^omehow special. It was initially re- 
ported to be CO,lI but later only A-type AFM order 
was found. □ However, Comparing the transport behav- 
ior of Fro.sSro.sMnOa with (La^Ndi_3)o.5Sro.5Mn03 of 
< z < 0.4, one finds similar behavior. Only the latter's 
resistivity levels off to a presumed metallic state slightly 
slowly after a jump at Tco or Tafm- So the bound- 



ary between the CO and A-type AFM states seems to be 
not so clear-cut: there may be a transition from the CO 
state to the A-type AFM state. For the most distorted 
class IV such as Pro.sCao.sMnOa and Ndo.sCao.sMnOa, 
no FM order appears .Q The PM phase changes directly 
into a CO state below Tco and then global AFM or- 
dering shows up at a lower temperature (~ 150K). The 
CO tendency in Class IV even extends to lower doping, 
although with|Ajpseudo-CE-type structure due to the ex- 
cess electrons .oil 
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FIG. 1. Experimental phase diagram of half-doped man- 
ganites. Lines are only a guide for eyes. Sources of the data 
are given in the text. The inset displays the elementary unit 
cell of Pnma and the numbering of the four Mn ions. 

The PM to FM transition is continuous. All the other 
transitions shown in Fig. yj are of first order with hys- 
teresis. The most striking feature of the CO statepis 
that it |ean be melt by rap external magncitic fieldtj, 
pressurea op-jelectric fieldtil or even by x-rayll3 or light 
irradiationgla into a FM state, indicating the competi- 
tion between them. The required magnetic field to melt 
the CO stateJo-Class IV is almost twice as large as that 
in Class III.E3^t£l In this paper, we shall concentrate on 
the peculiar phenomena associated with CO. The lattice 
structure of interest is orthorhombic with a space group 



Pnma (see Fig. |l| inset), which is well characterized and 
the most common in doped manganites. There are some 
scatters in the reported structures. This is understand- 
able because of the small distortion of the perovskitc 
structure, which is also sensitive to the preparation con- 
ditions. We develop below a Landau theory of the phase 
transitions in these systems and show that many struc- 
tural and thermodynamical behaviors are closely corre- 
lated within the framework. In particular, besides the 
magnetic, charge and structural/orbital pattern, primary 
features are a simple picture that unifies the three classes 
along with the nature of the transitions involved and the 
origin of the incommensurability of the orbital order, CE- 
type AFM order and the melting of CO state. 

We start with the magnetic transition. A FM phase 
transition is associated with a wave vector at k = or the 
r point. On the other hand, the CE-type AFM structure 
is described by k = (00^) and k = (^0^) for Mn ions at 
positions 1,2 and 3,4, respectively. In Lao.5Cao.5Mn03, 
the FM to AFM transition is found to be accompanying 
the incommensurate (IGj to (nearly) commensurate or- 
bital ordering transition.U In Class IV, the CE-type AFM 
order can appear separate from the charge and structural 
orders, which produce a natural configuration for the CE- 
type AFM structure (see below). Thus the AFM order 
is regarded as originating from the charge and orbital or- 
ders. For simplicity, we leave it towards the end of the 
paper. 

The PM to FM transition is described by an order pa- 
rameter M representing the average magnetization over 
the system. As pointed out in our previous work,t3 all the 
irreducible representations (IR's) at k = of Pnma are 
one dimensional. So we choose M as a scalar representing 
the compoaent that carries the IR r^ responsible for the 
transition.liZl Although it may couple to other magnetic 
configurations such as A-type or C-type AFM orders of 
the same IR, we can eliminate such modes if present and 
write the free energy simply as 



Fm = -aiM' 



\biM^ 



(1) 



with ai — aio{T — Ti), and aio and 6i depending only 
weakly on the temperature T. Eq. (w) describes a con- 
tinuous phase transition at T = Ti with M = yj —ai/bi. 

Next we move on to the CO phase transition. A promi- 
nent feature that indicates the existence of the CO state 
is the appearance of the superlattice diffraction spots 
characterized by a wave vector k = (^00) or the X point 
in the Pnma setting, which is adhered to throughout this 
paper. However, such spots associate more with lattice 
modulations than with the charge order that is character- 
ized by k = (100) or F. Direct charge and orbital orders 
are first detected by Xpsay resonant scattering techniques 
in doped La2SrMn04.ll3 So to describe the transition, we 
need two order parameters. 

On the one hand, the charge order is characterized 
by an order parameter C = ^i + ^2 ~ ^3 ^ ^4 with S,i 



the occupancy possibility of site i (Fig. |l] inset). C is 
maximized if ^1 = ^2 — ^^3 = ^^i, and so a nonzero 
C produces the observed CO pattern. By considering 
the permutations of the four sites under the symmetry 
operations of the Pnma group, it can be shown that C 
transforms also as the IR r^ of Pnma at F. Accordingly, 
the free energy of the CO transition is given by 



Fc = ^a2C^ 



:b,C' 



(2) 



where 02 — a2o(T — T2), and a2o and 62 are constants. 
Below T2 charge order appears^apd the symmetry of the 
structure is lowered to P2i/mS3 

On the other hand, the structural modulation at 
k — (^00) is described by one of the two IR's of 
the wave vector —namely, Xi and X2, which are both 
two dimensional. EZI Accordingly, the structural transition 
is characterized by a two-dimensional order parameter 
(^17^2)- The physical meaning of the order parameter 
may be the displacement of the Mn^+Oe octahedra as 
modeled by RadaelU and coworkers to account for the 
diffraction patternsB It can be shown that 772 rnay rep- 
resent an arbitrary linear combination of the x and z 
components and 771 the y component of such displace- 
ments. Such a displacement pattern is consistent re- 
spectively with an orbital configuration of d^^2_j.-2 and 
dx2-y2 , which, when propagate half a period to X, switch 
to d^z^-r^ and dy2_z2, respectively. It may be possible 
to choose alternatively an orbital basis such as d^^2_j.2 
and dy2_z2 and then represent the order parameter as 
the long-range order of a certain orbital which assumes a 
certain angle in the orbital space. We just note in passing 
that a certain displacement pattern corresponds to some 
orbital order. 

Notice that the frequently observed P2i/m symmetry 
can only arise from the IR ATiEj. Therefore, the free 
energy for this structural transition is 



F, 



1 



1, 



1 



where again 03 — 030 (T — r3), and 039, &3, d, k and a 
are constants. 

A peculiar feature of Eq. (0) is the appearance of the 
Lifshitz invarianLiilie k term), which frequently leads to 
IC modulations.E30 Many characteristic features of an 
IC transition has been observfid—in the CO and struc- 
tural transition in manganites .aE3 Therefore, the IC na- 
ture of the structural/orbital (but not charge) order has 
its origin in the Lifshitz invariant of the X point. Nev- 
ertheless, we shall neglect this IC feature of the struc- 
tural modulation below for simplicity and focus on its 
interplay with the charge and magnetic orders. This is 
partly justified by the fact that commensurate structure 
is also frequently observed in the same experiments that 
display the reversefl In this case, Eq. (||) then exhibits 



two possible phases belo w Tj,. One ha s only one nonzero 
component ecuial to ±^J —a^/{h-i + d) and so its symme- 
try is P2i/niB if 63 + d > and 63 < 0. The other 
satisfies 771 — ±j]2 = ±^—03/(63 + 2d) and belongs to 
Pra symmetrytj when 63 + 2(i > and 63 > 0. 

Coupling of the charge to the structural degrees of free- 
dom can be readily found by noting that Tq\ — rj^ trans- 
forms as the same IR as C, and so the simplest coupling 
between them is 



Fcv = 9CvC{vf -~vi), 



(4) 



where gcri is a measure of the coupling. Assume gcr/ and 
C are positive without loss of generality. It is transparent 
then that the coupling will favor the ordering of 772 once 
the charge is ordered, since the transition point for 772 is 
now elevated to 03 — 2gcr]C, while that for 771 lowered to 
«3 + '^9CriC- This explains the reason for the often ob- 
served P2i/m symmetry.cl Possibility for ordering of both 
?7i and 772 still exists, which may account for the absence 
of the 2i screw axis in Smo.5Cao.5Mn03, which should 
then be of Pm synametry, though it was preferred to be 
P2mm or Pmmrn.c3 This scenario is confirmed from the 
phase diagram illustrated in the inset of Fig. 0. It is seen 
that as 63 increases, the boundary of the Pm phase moves 
to the right, reducing the region of the P2i/m phase. If 
&3 < 0, only the P77T- phase exists, while for d < 0, the 
Pm phase extends far to the right. Note that a large 63 
or small d means a large "lock-in" term [the 63 term m 
Eq. (P)] that tends to suppress the incommensurability.Ell 
This seems to be consistent with the observation that 
more distorted systems such as Sm and Gd, tend to be 
more stoichiometric and so commensurate.B^I We remark 
that the transition from P2i/m to Pm phase may accom- 
pany or be hidden by the IC to commensurate transition 
and PM to CE-type AFM transition. 

Since AI changes sign by time reversal, the only possi- 
ble couplings of the magnetic to the charge and structural 
transition are bi-quadratic, i.e.. 



F, 



MCri 



^9Mr,M^iv! 



V2) 



-QMCC^M^, 



(5) 



where both coupling coefficients are positive due to the 
competing orders. 

Equations (|l[]|) constitute our theory of the magnetic 
and charge and structural or orbital transitions. Instead 
of going into detailed estimations of the various coef- 
ficients in the model, we are content here with global 
features that are believed to be relevant to the parame- 
ter regime of real materials. To this end, we study of a 
simplified version of the theory, in which we have taken 
771 = 0, i.e., disregarded the possibility of the Pm symme- 
try, relabeled b^ + d as 63, and neglected the bi-quadratic 
coupling between the charge and magnetism, since their 
coupling to the structural order results in a lower-order 
CAf ^-type coupling. This implies that long-range charge 
order always accompanies with structural order. 
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FIG. 2. Theoretical phase diagram of the coupled mag- 
netic, charge and structural transitions. The parameters are 
aio = a20 = 0,30 = &i = 62 = 63 = 1, Ti = 3, T2 = 2.5, T3 = 2 
and QMr] = 1.2 (solid lines). The dashed line corresponds to 
(7M17 = 1.5, while the dotted lines to Ti = 2.5. Inset: phase 
diagram of the charge-structural transition. The dashed lines 
represent bs — 2 while the dotted ones d = —0.3. 
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FIG. 3. Order parameter vs temperature T . Thick lines 
correspond to gcr, = 1.3 while thin and mediate lines gcr, ~ 1 
and gcri = 0.3, respectively. The inset shows the correspond- 
ing free energy from bottom. The upper curve is for the pure 
FM phase [Eq.(|l)j. Other parameters are the same as the 
solid lines in Fig. 3. Overlaps in gcri ~ 0.3 imply metastable 
region. 

Figure ^ displays a generic phase diagram for several 
sets of the parameters in arbitrary units. It shows the 
strong coupling case gj^j > bibs between magnetism and 
lattice so that no mix magnetic and charge orders ap- 
pears. It can be seen that as gcrj increases, the sequence 
of phase transitions changes from a PM -^ FM to a PM 
^ FM ^ CO then to a PM -^ CO phase directly as T 
is lowered. The real situations as in Fig. Ill are of course 
not simply only a variation of gcrj, but Fig. exhibits in 
a simple way the relevance of the theory. Figure ^ shows 
the variation of the order parameters for three different 
gcri- It is seen that the transition to CO phase is dis- 



continuous. The variations of C and VrSire similar, but 
their couphng is not necessarily linear.O Note that the 
CO and the structural transitions are so strongly coupled 
that they take place at a single transition temperature. 
Although the CO can still appear preceding the struc- 
tural order when gcri < 1, the reverse is not true, namely, 
77 = if C = 0. Therefore the structural or orbital order 
is driven by the charge order. The inset plots the cor- 
responding free energy vs temperature, showing that the 
larger the gc?)7 the lower the free energy, and so the bigger 
the magnetic field that lowers the upper curve [Eq. (|^)] 
by HM to melt the CO state, being in agreement with 
experiments. This indicates that the magnetic field acts 
in melting the CO state more than the AFM state. An 
evidence for this is that above the CE-type AFM order- 
ing temperature, no AFM for the field to melt. This also 
justifies the separate treatment of the AFM order. 

Finally we discuss briefly the transition to the CE- 
type AFM order. Noting that in Lao.5Cao.5Mn03 with 
the standard CE-type state, the magnetic moments lie in 
the a-c planen, we may choose the x and z components 
of the AFM vectors Li = /Xj^ — /Xj and L2 = /I3 — /i.4 
as order parameters, since they transform as the two 
components of the two-dimensional IR's r^ and r^ com- 
bining with its complex conjugate r'^ associatecLrespec- 
tively with the wave vectors (00 5) and (^0^),EZI where 
/Xi is the magnetic moment of ion i. The two IRIs. 
rather than one makes the transition discontinuous .Ell 
Both IR's give a magnetic symmetry of Pb2i/m and a 
structural one of P2i/m which is identical with the CO 
and structural transitions. Ej The compatibility of the 
orbital and magnetic patterns implies an enhancement 
of both transitions and can be described by a coupling 
(j]l + ri2){LiaLif3 + L2aL2i3) with a negative coupling con- 
stant, where a and /3 denote x or z. As a result, ordering 
of one kind of the orders enhances the other, leading, for 
instance, to the accompanying of the AEM order with 
the IC to commensurate orbital ordeiaE3 as the AFM 
order promotes the contribution of the lock-in term, and 
to the increasing of resonant x-ray scattering-intensity of 
charge and orbital orders as AFM orderingo 

In conclusion, we have developed a Landau theory for 
the coupled phase transitions in half-doped manganites 
through the identification of the order parameters for the 
FM, CE-type AFM, CO and structural or orbital orders. 
The theory provides a unified picture for the scenario 
of the phase transitions and their nature with respect 
to the variation of the tolerance factor of the manganites 
via the symmetry-adapted coupling among the degrees of 
freedom. Many peculiar phenomena of half-doped man- 
ganites result from the interplay between the FM or A- 
type AFM and CO states, the CE-type AFM sets in only 
as a secondary factor. So an applied magnetic field pri- 
marily melts the CO state. The theory also accounts 
for the origin of the IC nature of the orbital order and 
its subsequently accompanying AFM order. As a phe- 
nomenological theory, it can make direct contact with the 
experimental results especially the symmetry of the in- 



volved structures that sensitively infiuence the transport 
behavior in manganites. Experimental clarifications are 
desirable of the symmetry of the CO state and its rela- 
tion to the oxygen stoichiometry and commensurablity, 
possible structural transition from P2i/m to Pm, and 
the relation between A- type AFM and CO states. 
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